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mSITU  CHARACTERIZATION  OF  SATURATED  SANDS  AND  SILTS  FOR  THE 
PREDICTION  OF  DYNAMIC  SHEAR  MODULUS  AND  SHEAR  WAVE  VELOCITY 


1.  INTRODUCTION 

The  majority  of  the  numerous  analytical  methods  presently  available  for  assessing 
the  response  of  soil  deposits  or  soil-structure  systems  to  earthquakes,  explosives  or 
machine  loadings  require  an  assessment  of  the  shear  modulus  G  at  shear  strain 

fDdX 

amplitudes  less  than  .001  percent.  Field  values  of  are  used  directly  to  determine 

soil  stiffness  for  low  amplitude  vibration  studies  (Rlchart  et.al.,  1970).  For  high 

amplitude  studies,  field  values  of  are  used  as  reference  values  (Seed  et.al.,  1970). 

rriax 

Geophysical  methods  utilizing  cross  hole  or  down  hole  technique  are  often  used  for 
critical  soil  analyses.  However,  in  certain  situations  empirical  equations  or  laboratory 
testing  methods  are  used  to  determine  G  . 

illoX 

Empirical  equations  used  for  sarxls  are: 


^max 

1230  -  (o  )^ 

i+e  o 

Hardin  and  Black  (196S,1969) 

(1) 

^max 

‘OOO  ■'amax 

Seed  and  Idriss  (1970) 

(2) 

G  = 

max 

1220  N°** 

Ohsaki-Iwasaki  (1973) 

(3) 

Where  e  is  the 

void  ratio,  a  is  the  mean 
o 

effective  confining  pressure  in  psi  for 

the 

Hardin-Black  equation,  and  in  psf  for  the  Seed-Idriss  equation,  K2max 


is  an  empirical 


factor  which  varies  according  to  the  relative  density,  N  is  the  standard  penetration 

2 

value  and  G _ ,  in  the  Ohsaki-Iwasaki  equation,  is  in  t/m  . 

max 

Equation  (3)  suffers  limitation  associated  with  an  erratic  procedure.  In  particular, 
the  values  of  N  can  vary  as  much  as  fifty  percent  and  hence  it  is  not  considered 
further  in  this  study. 

Equations  (1)  and  (2)  give  similar  results;  Seed's  and  Idriss*  procedure  for  the 
verification  of  equation  (2),  used  in  part  the  data  developed  by  Hardin  and  Black. 
Equation  (2)  will  be  considered  for  further  laboratory  correlation. 


Investigations  have  shown  that  modulus  values  for  sands  are  strongly  influenced 
by  the  confining  pressure,  the  strain  amplitude,  the  void  ratio,  and  the  angularity  of 
the  particies  at  low  confining  pressures;  they  2ure  however  not  significantly  affected 
by  grain  size  characteristics.  Very  limited  data  is  available  to  indicate  the  influence 
of  anisotropy  on  Thus  it  can  be  argued  that  the  parameter  K2max  equation 

(2)  is  a  function  of  porosity,  shape,  and  orientation  of  particles  and  that  G  is  a 
direction  dependent  property.  *^2max  depends  on  cementation  of  particles  as 
pointed  out  by  Seed  and  Id’iss  (1970).  Thus,  to  completely  characterize  a  sand  taking 
into  consideration  porosity,  shape,  orientation  of  particles,  anisotropy  and  cementation, 
a  methodology  is  needed  so  that  values,  obtained  from  in  situ  shear  wave 

velocities  measurements,  can  be  correlated  with  the  structural  characteristics  to  predict 
shear  wave  velocity.  It  is  the  purpose  of  this  study  to  use  the  recently  developed 
methodology  of  characterizing  sands  by  electrical  methods  to  develop  such  a  correlation. 

2.  NON  DESTRUCTIVE  METHOD  TO  CHARACTERIZE  GRAIN  AND  AGGREGATE 
PROPERTIES  OF  SANDS  AND  SILT 

It  is  generally  accepted  that  the  derived  properties  of  sand  deposits  are  governed 
by  the  grain  and  aggregate  characteristics  of  the  particles.  It  has  been  shown 
experimentally  by  numerous  studies  carried  out  by  Archie  (1942),  Jackson  (1975),  Wyllie 
and  Gregory  (1953),  Arulanandan  (1975),  Aruianandan  and  Kutter  (1978),  Kutter  (1978), 
and  Arulrholi  (1980),  that  a  nondimensional  electrical  parameter  called  the  Formation 
Factor  is  dependent  upon  particles'  shape,  long  axis  orientation  contact  orientation 
and  size  distribution  and  also  void  ratio,  cementation,  degree  of  saturation  and 
anisotropy.  The  formation  factor,  F,  is  defined  as  the  ratio  of  the  conductivity  of 
the  electrolyte,  o^,  to  the  conductivity  of  the  saturated  sand-solution  mixture,  o^. 
The  dependence  of  formation  factor  on  void  ratio,  particles'  shape  and  long  axis 
orientaton  has  been  theoretically  shown  by  Arulanandan  and  Daf alias  (1979)  and  Daf alias 
and  Arulanandan  (1978).  It  has  also  been  shown  from  thin  section  studies  and  formation 
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factor  measurements  in  different  directions  that  the  structure  of  sand  maiss  in  most 
cases  is  transversely  isotropic,  with  the  vertical  axis  being  the  axis  of  symmetry,  both 
in  natural  deposits  and  in  samples  prepared  in  the  laboratory.  However,  in  certain 
situations,  the  manner  in  whidi  the  sand  deposit  is  formed  may  cause  the  sample  to 
be  orthotropic.  But  to  a  close  approximation,  it  is  reasonable  to  consider  the  structure 
of  the  sand  mass  to  be  transversely  isotropic. 

The  formation  factor,  when  suitably  interpreted,  can  be  used  to  quantify  and 
predict  the  porosity,  shape  and  anisotropy  of  sand  deposits.  For  example,  an  average 
formation  factor  F  given  by 

F  =  (Fy  +  2F^)/3 .  (<») 

with  Fy  and  F^  being  the  vertical  and  horizontal  formation  factors  respectively,  has 
been  shown,  both  by  theory  and  by  laboratory  measurements,  to  possess  the  properties 
of  the  first  invariant  of  a  second  order  tensor.  A  tensor  is  defined  as  a  quantity 
having  physical  significance  whose  components  with  respect  to  a  coordinate  system 
transform  according  to  a  specific  law  upon  change  of  the  coordinate  system.  A  second 
order  tensor  has  nine  components.  Tensor  transformations  are  linear  and  homogeneous 
and  hence  if  tensor  equations  are  valid  in  one  coordinate  system,  they  are  valid  in 
any  other  coordinate  system.  This  invariance  of  tensor  equations  under  a  coordinate 
tranformation  is  one  of  the  principal  reasons  for  the  usefulness  of  tensor  methods  in 
applied  mechanics.  It  has  been  shown  that  F  is  independent  of  anisotropy  caused  by 
preferred  particles'  orientation  and  is  a  direct  measure  of  porosity,  n,  for  a  given  sand. 
To  illustrate  this,  horizontal  and  vertical  formation  factor  measurements  on  lucite  balls 
and  Monterey  '0/30*  sand  samples,  prepared  by  three  different  methods,  were  made  in 
laboratory.  Monterey  *0/30*  sand  is  uniformly  graded  and  lucite  balls  are  1/8  inch 
diameter  spheres.  Two  six  inch  cubical  cells,  one  with  two  six  inch  squares  platinum 
coated  copper  electrodes  fixed  on  two  opposing  vertical  faces  and  the  other  with 
electrodes  on  the  top  and  bottom  faces  have  been  used  for  making  horizontal  and 


vertical  electrical  resistance  measurements  respectively.  The  formation  factor  values 
were  calculated  using  the  following  equati<Mis; 


F  = 


a  - 

s  R  A* 


where  R  is  the  resistance,  L  is  the  distance  between  the  electrodes  and  A'  is  the  area 
of  the  electrodes.  The  resistance  measuremmts  were  made  to  half  a  present  accuracy. 
F-n  relationships  for  all  three  methods  were  seen  to  be  identical  and  are  shown  in 
Figs.  1  and  2  for  lucite  balls  and  Monterey  '0/30'  sand  respectively.  Five  more  uniform 
sands  and  a  silty  sand  from  the  Revelestoke  dam  site  were  also  used  in  the  study  and 
unique  relationship  between  F  and  n  was  found  to  exist  for  each  sand.  The  F-n  curves 
for  all  the  sands  and  lucite  balls  are  shown  in  a  Log-log  plot  in  Fig.  3.  It  would 
appear  that  the  average  formation  factor  for  any  type  of  sand  describes  the  state  of 
packing  of  sand  and  is  independent  of  the  long  axis  orientation  of  the  particles. 

The  unique  relationship  between  F  and  n  can  be  used  to  predict  in  situ  porosity 
of  uncemented  sands.  The  accuracy  of  the  use  of  this  method  for  in  situ  prediction 
of  the  porosity  of  cemented  sarxis  needs  to  be  investigated. 

2.1  Shape  Factor  and  Anisotropy  Index 

An  integration  technique  proposed  by  Bruggeman  (1935)  was  used  by  Dafalias 
and  Arulanandan  (1975,1978)  to  derive  an  expression  for  avereige  formation  facotr,  F, 
as  a  function  of  porosity,  n,  and  average  shape  factor,  f,  as 


F  =  n" 


The  average  shape  factor  f  b  the  negative  slope  of  the  log  F-Log  n  plot.  It  is  the 
first  invariant  of  the  second  order  shape  factor  tensor  f  and  it  relates  the  electric 


fields  inside  and  outside  the  sand  particles.  It  has  been  shown  both  theoretically  and 
experimentally  that  the  shape  factor  is  directional  dependent  and  depends  on  porosity, 
gradation  and  particles'  shape  and  orientation  Arulmoli  (19S0),  Daf alias  et  al  (1979), 
Kutter  (1978).  Since  the  average  formation  factor  is  independent  of  orientation  of 
particles,  the  average  shape  factor,  for  a  given  sand,  is  expected  to  be  a  function  of 
porosity  and  the  shape  of  particles.  A  mean  value  of  the  average  shape  factor,  f^ 
defined  as 


mean’ 


^mean  "  2  ^*max  ^min^ 


(8) 


where  f^^„  and  f_,.  are  the  extreme  values  of  f  at  extreme  porosities,  is  used  in 
max  min  ^  ’ 

the  correlations  proposed  in  this  paper.  The  maximum  percentage  difference  between 

and  ?  is  about  five  percent,  Arulmoli  (1982). 
max  min  ^ 

An  anisotropy  index  A  was  introduced  by  Arulanandan  and  Kutter  (1978)  as 


A  =  (Fy/Fj^) 


1/2 


(9) 


Anisotropy  of  sand  structure  is  due  to  the  orientation  of  individual  particles  and  the 
contact  orientation.  Formation  factor  measurements  made  on  samples  consisting  of 
spherical  particles  showed  anisotropy  for  certain  structural  arrangements.  Even  though 
spherical  particles  do  not  have  any  preferred  orientation,  the  contact  orientation  was 


significant  enough  to  produce  anisotropy  in  structure.  Formation  factor  measurements 
are  thus,  sensitive  to  contact  orientation  of  particles,  in  addition  to  the  orientation 
of  individual  particles.  This  has  been  elaborated  theoretically  in  detail  by  Dafalias 
and  Arulanandan  (1979).  Dependence  of  formation  factor  on  cyclic  stress  history  has 
also  been  shown  by  Arulmoli  (1980).  If  the  samples  of  sands  are  consolidated  under 
different  anisotropic  stresses,  it  is  reasonable  to  assume  that  the  formation  factors 
will  be  different  for  different  stress  conditions.  This  difference  is  expected  to  be  due 
to  the  particle  and  contact  orientations  being  different  from  different  stress  conditions. 


2J2  Overall  bidices  of  Soil  Characteristics 


It  is  seen  from  Fig.  4  that  the  derived  porosities  of  a  sand  deposit  are  governed 
1^  the  grain  and  aggregate  characteristics  of  the  particles.  It  is  clear  from  the 
foregoing  paragraphs  that  the  formation  factor  depends  on  these  characteristics  and 
hence  electrical  parameters  obtained  usng  formation  factors  of  sands  in  different 
directions  may  be  used  to  correlate  with  soil  properties  such  as  liquefaction  potential, 
friction  angle,  permeability  and  compressibility^  The  parameters  F,  A  may 

be  combined  to  obtain  such  correlations.  These  empirical  correlations  would  be 
extremely  useful  in  evaluating  the  performance  of  sites  which  contain  sand  deposits. 

2.3  Electrical  Probe 

Field  values  of  formation  factors  were  obtained  using  an  electrical  probe. 
Geoelectronics  Model  GE-100,  shown  in  Figure  5  (Arulanandan  (1977)).  The  probe 
consists  of  three  main  parts: 

1.  The  main  body  is  a  4  feet  long  hollow  steel  tube  with  a  3  inch  outside  diameter 
and  a  wall  thickness  of  1/8  inch.  An  electrical  bridge  is  housed  inside  the  main 
body  together  with  a  12  volt  gear  pump  to  pump  solution  into  a  cell  for  solution 
conductivity  measurements.  The  bridge  can  resolve  resistance  to  an  accuracy  of 
about  one  percent. 

2.  The  probe  tip  is  a  12  inch  long  steel  tube  with  a  3  inch  outside  diameter  and  1/16 
of  an  inch  wall  thickness.  The  probe  tip  carries  the  electrodes  for  sample 
measurements  in  two  different  directions.  Figure  6.  A  tiny  open  cell  with  electrodes 
is  also  located  along  the  outside  wall  of  the  probe  tip  and  it  is  covered  with  a 
porous  stone.  The  probe  tip  can  be  unscrewed  from  the  main  body  and  is  replaceable. 

3.  The  control  unit  consists  of  a  microprocessor  for  transmitting  the  electrical  signal 
and  receiving  the  response  of  the  material  within  the  probe  tip  in  the  form 
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of  resistances  and  capacitances.  The  controi  unit  is  connected  to  the  electronics 
in  the  nn2un  body  through  a  150  foot  long  stiff  cable.  The  whole  probe  assembly 
can  be  operated  with  a  12  volt  d.c.  power  source. 


3.  LABORATORY  CORRELATION  BETWEEN  ELECTRICAL  PARAMETERS  AND 
SHEAR  MODULUS 

One  way  of  predicting  maximum  shear  modulus,  G  ,  is  by  measuring  the  in 

IKlaX 

situ  shear  wave  velocity,  V^,  and  by  using  the  equation 


G  =  p  V 
max  ^  s 


where  p  is  the  mass  density  of  the  deposit  at  the  depth  of  measurement.  Investigations 
have  shown  that  the  maximum  shear  modulus  values  for  sands  are  strongly  influenced 
by  the  confining  pressure  and  the  void  ratio.  A  relationship  between  the  shear  modulus, 
G,  in  psf,  and  mean  effective  confining  pressure,  o^,  in  psf,  was  given  by  Seed  and 


Idriss  (1970)  as  (2): 


G  =  1000  K,(o  . (11) 

i  o 


where  the  parameter  K2  depends  on  the  void  ratio,  strain  amplitude,  geological  age 
of  the  sand  mass  and  in  situ  stresses.  Thus,  the  maximum  shear  modulus  at  very  low 
shear  strain  amplitudes  is  related  to  through  the  maximum  value  of  K2, 

and  is  given  by  equation  (2) 


=  1000  3  . <2) 

^max  ® 

In  the  equation  (2),  K2max  largely  on  void  ratio  and  also  on  the  age 

of  the  deposit.  Field  and  laboratory  shear  wave  velocity  measurements  were  made 
using  cross  hole  and  down  hole  seismic  measurements  and  the  electrical  measurements 
were  made  in  the  field  and  in  the  laboratory  using  the  electrical  probe.  Shear  wave 
velocity  measurements  in  the  field,  were  made  at  three  sites:  El  Centro  Site  (Ertec 
report  (1982)),  Lawsons's  Landing  Site  at  San  Francisco  (Kleinfelder  and  Associates 


(19S1)),  and  Tientsin  Site  in  Peopie's  Republic  if  China  (Arulanandan  (19S2)),  and 
laboratory  measurements  made  on  Ottawa  Cl 09  and  soil  from  Lawson's  Landing  site 
were  used  to  estabiish  a  correlation  between  obtained  from  equation  (2)  above 


—  —  1/2 

and  an  electrical  parameter.  The  electrical  parameter  F/(A«f^)  was  found  suitable 


to  correlate  with  K. 


‘2max‘ 

The  field  results  are  the  laboratory  results  of  shear  wave  velocities  and  the 


electricai  properties  are  given  in  Tabie  1,  2,  3  and  4.  The  gradation  characteristics 
of  the  sands  and  the  iaboratory  vaiues  of  electrical  properties  are  given  in  Figs.  7,  S, 
and  9.  The  correlation  between  and  F/(A*f^)  is  given  in  Fig.  10. 

In  situ  measurements  of  the  formation  factors  in  the  horizontai  and  vertical 
directions  enabie  us  to  determine  F,  and  A  and  laboratory  determination  of  F-n 
relationship  would  give  us  hef'C®  *^2max  ^  predicted  using  the  reiationship 


shown  in  Fig.  10.  The  value  of  can  be  <'alculated  from  equation  (2)  for  the 

MlClX 


depth  of  the  deposit  under  consideration. 

A  typical  example  of  the  comparison  between  the  measured  and  predicted  shear 


wave  veiocities  with  depth  using  the  above  procedure  is  shown  in  Fig.  11. 
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